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ABSTRACT 


In atmospheric experimentation three ground-based detectors suffice 
for determining wind vectors that could be present in a certain space 
volume. While this is true in general, the geometry simplifies a great 
deal, when these vectors can be taken as horizontal. Since to do so will 
often be warranted, the pertinent method for a proper arrangement of 
detector positions and attitudes has been developed in detail. With a 
suitable layout winds reasonably constant near a selected height, if 
blowing into and out of a quarter of the rose, can be calculated with 
sufficient accuracy and confidence. The error analysis assumes that 
the observational time uncertainty is at most ±0.1 sec. 

To observe wind vectors in the remaining directions, or at several 
heights simultaneously, the use of more detectors (or of the available 
multiple detectors) is called for. This is subject matter for future 
investigation, as is the more intricate problem posed by non-horizontal 
winds . 
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TECHNICAL MEMORANDUM X- 53754 


WIND VECTOR CALCULATION USING CROSSED-BEAM DATA 
AND DETECTOR ARRANGEMENT FOR MEASURING HORIZONTAL WINDS 

SUMMARY 


The report's first part (sections I through V) contains the geo- 
metric foundation and subsequent development of a general method for 
measuring atmospheric winds by exploiting the ideas of crossed-beam 
experimentation. It shows that (and how) with three single-beam 
detectors at ground level one should be able to monitor a variety of 
unknown winds that could be blowing within a limited space volume, pro- 
vided any such wind, once present, remains sufficiently stable during 
observation. 

A second part (sections VI through VIII) applies -- and further 
develops -- the method in the practically meaningful and geometrically 
least complex case where the winds can be taken as horizontal. A cri- 
terion for this to be justifiable is given. Flow conditions, wind 
identification, error analysis, and the desire to keep the detector 
arrangement compact impose a number of restraints that are reflected in 
the experimental layout finally suggested. These restraints are more 
acutely felt the higher up one proposes to measure. Equations for 
actually computing the two velocity components are supplied. 

A concluding section delineates problems yet to be solved. 


I. INTRODUCTION 


As it is well known, wind determination by crossed beams presents 
a number of problems. By way of examples we may list: record evaluation 

complications through long observation times and non-s tat ionary statistics, 
interacting intensity fluctuations through the scattering of light from 
different sources (sky and sun), the disturbing influence of cumuli and 
other clouds, unwanted features in absorption physics. Most of them are 
or have been under vigorous attack elsewhere. 

The report on hand is concerned with the specific difficulties 
encountered when the net results of the experimentation proper, i.e., 
the characteristic delay times, are to be used for calculating the winds 
they are associated with. 



To achieve this end, resources of spatial geometry are assembled in 
a preparatory section. Still there are quite a few problems left that 
must be overcome. In the first place, covariance peaks at characteristic 
times will be registered whenever eddy trains travel on paths that happen 
to connect beams; one must sort out those that belong to one and the 
same wind. In the second, one should be able to say in what space region 
the latter was blowing. Thirdly, one cannot be satisfied to detect just 
this one wind on the chance it might be present, but wishes to monitor 
an array, as large as possible, of unknown wind vectors near a predeter- 
mined location. Fourth, observa tional error margins must be prevented 
from getting blown up to unacceptable dimensions. For this reason travel 
paths too short cannot be admitted. Too long ones are prohibited lest 
eddies lose identity. Finally, one would not like, on practical grounds, 
to have to place the detectors too far apart. 

These warnings and requirements were heeded in working out the body 
of the paper. A rather definite set of detector specifications for 
observing horizontal winds is the outcome that leaves to free choice but 
two parameters out of ten. 


II. GEOMETRIC PRELIMINARIES 


Geometric analysis provides a number of results useful for the 
argument to follow and for the realization in general of spatial rela- 
tionships associated with the crossed-beam method. 

The two straight lines, a and b, established by the points P P 2 
and the direction vectors*, 

a = a]± + 

j3 = Pii + Pal + fell* 


Vectors are designated by underlined symbols; e.g., i^, k are the 
unit vectors in axis directions. 



have the equations 


i.e., that 
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This excludes parallelism as well. Nevertheless, the lines do define two 
parallel planes in which they are contained. The direction of that pair's 
common normal is given by the unit vector (called the "binormal" direction 
here*) 


n = — 7 
— sm u) 


ab 


[a X _g] , 


( 2 ) 


u> a b being the angle made by the vectors CL and when crossed: 


cos w ab = aiPl + a + <*3fo- 


Introduce now a further direction 


m = m x i_ 4- m^i + m 3 k ( | m = 1 ) 


This agrees with the usage of the word in the theory of spatial curves, 
provided that the two planes are interpreted as the lines' tangential 
planes . 
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not parallel to those planes: 



m i 

m 2 

m 3 

(map) = 

a x 

a 2 

a 3 


Pi 

P2 

P? 

where m = 

= n 

evidently 


0 . 


Then there exist two -- and only two -- points, P-j_ and P 2 , on the 

lines a and b such that P* P* is parallel to m. Clearly, among the 
directions m only those are admitted as are capable of connecting a with 
b in that order (excluding all the reverse orientations). The position 
vectors, and rjp, of P* and P 2 are given* by 
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= r 


i + (map) 
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— — + (map) 
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z 2 - z x 

<*1 

a 2 

a 3 
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m 2 
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> , ( 3 ) 


where r x and r 2 define P x and P 2 . It is clear and can be shown mathe- 
matically* that the choice of these two points on a and b cannot affect 
the vectors r* and r 2 . If the length P* P 2 is denoted by we may write 


p* p* 


• • * X") m 

r 2 " £i - R ab S* 


Figure 1 illustrates this formula. 


Vc 

See: Heybey, W. H. , "On the Wind Component Measurable by Crossed-Beam 

Arrangement," IN-AERO-7-67 , October 12, 1967. 
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Figure 1. The vectors r*, r*, and m 


The projection onto the binomial of any straight segment connecting 
a and b is equal to the binormal distance, 5, itself. This is readily 
recognized by figure 2, which depicts the simplest geometric arrangement 
possible without loss of generality. Of the two parallel planes estab- 
lished by the lines a and b one is taken as the (z,x) -plane and the 
other as parallel to it. The y-axis is defined by the binormal. The 
projecting lines through the points Pi or P* on the one hand, and through 
or P 2 on the other are the lines a and b themselves, as they are normal 
to the binormal. 


•JU JU 

Since thus the projections of P" P 2 and Pj P 2 are equal: 
(r'| - r*) • n = (r 2 - r x ) • n. 
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Figure 2. General two-beam geometry. 


we find that 
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or that 



This is the shortest expression for the distance P* p*. That the 
quotient at right is always positive can be verified in a convenient 
manner if (as one is permitted to do) the defining points P l5 P 2 are 
replaced by the binormal terminals A, B with coordinates 


x i = yi = zi = 0 

x 2 = 0, y 2 = 6 > 0, z 2 = 0. 


Then, since a 2 = p 2 = 0, 


A ab 


0 & 0 
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(map) = 
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Because the vector m points from to P 2 , its y-component is positive: 
m 2 > 0. It follows that 


A 


‘ ab = — > 0. 


(map) m 2 


III. TWO- BEAM EXPERIMENTATION 


Broadly speaking, wind determination by crossed beams assumes that 
the bulk speed of eddy propagation can be taken as a sufficiently close 
approximation to the main wind speed. In the simplest atmospheric 
experimentation two detectors receive natural light out of a conical 
volume extending from their objectives to great distances. For the 
present purpose these "beams" are replaced by the axes of the cones. 
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Likewise, turbulent eddies will be idealized to have infinitesimal size. 

The time, an eddy needs to cover the distance between two points 

Pi and P* is determined experimentally. It is assumed as equal to that 
particular delay time that produces maximum covariance between the intens- 
ity curves obtained from the two beams, and is positive or negative depend- 
ing on whether, in order to locate the maximum, the b-record has to be 
delayed relative to the a-record or vice versa. If it is negative, the 
velocity 


V = 


R 


* 


ab 

T * iii 
ab 


m 


(5) 


goes from P* to P*, i.e, it is in the reverse direction of m. 

Wind variations pending observation are to be reckoned with. Fluc- 
tuations in strength tend to blur the covariance maximum, but would not 
blot it out completely except perhaps when excessive or unidirectional. 
Variation of direction appears more serious. The eddies leaving P'£ 

(or P*) over a period of time may often miss the second beam*; those that 
do arrive at it may not be numerous enough for a clear definition of the 
maximum. Although the beams a and b are immersed in a three-dimensional 
wind field, the directions of the motion could be inappropriate at every 
point of the a- (or b-) beam, meaning that no single eddy path will cross 
the second beam; compare Figure 3 (where however two connecting paths 
are indicated). At any rate, the number of covariance maxima will be 



Figure 3. Curved eddy "sheet" ABCD downstream from beam a. 


*In physical reality this danger abates somewhat because of the finite 
beam width and finite eddy size. 

8 



limited and may even be zero, for instance, if the wind is blowing in 
a direction normal to the beams 1 binormal. This observation is a first 
indication that more than just two detectors will be needed for meas- 
uring winds. 

Counting on the finite beam width and eddy size, we infer from 
Figure 3 that it will be advantageous to have the beams pass each other 
at a relatively close distance" to insure eddy interception and an 
approximate wind vector finding. However, direct intersection would 
cause maximum correlation at the common point, irrespective of wind 
direction (which thus remains unknown) and leaving us with an indeter- 
minate expression for the wind strength, V, the distance zero being 
covered in zero time."* A finite binormal distance is required which 
moreover ought not to be chosen too small in order to escape short 
travels with consequent large effects of observational inaccuracies. 

On the other hand, it must not be allowed to be too large either, lest 
eddy change or decay on the longer paths impair or destroy the needed 
correlation. 

In an ideal atmosphere (considered first) the wind is constant 
within sight, its direction being given by the unit vector 


v = v x i *f + v 3 k. 


( 6 ) 


If we knew it, we should put 


m = ± v 


(7) 


and, by the system (3) could uniquely determine the position vectors 
r* and v% (the double sign cancels). The covariance curve would 
exhibit” single maximum at the time, needed for an eddy to move 

from P* to p| (or conversely). The strength of the wind would emerge as 



A later reasoning will throw support to beam nearness from quite a 
different angle. 

If, in a rare instance, the wind vector happened to be parallel to 
the plane of the intersecting beams, the situation would be quite 
as confused. 
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This result applies in wind tunnel experimentation when the convective 
motion can be taken as parallel to the tunnel axis. In the atmosphere a 
reasonable surmise regarding v will not in general be at hand. However, 
the component of the wind vector 


V = iVi + jV 2 + kV 3 


( 8 ) 


in direction of the binormal is still obtainable. With the aid of 
expressions (5), (7), (4), and (2) 


„ _ . Jib . , = + ^ab S- • n -ab (vaB) 1 

n ” x* b ' x* b (± vaB) T^ b (vccb) sin 


u) 


ab 


or 


A ab 1 

V = “4 : . 

n t - sin oj , 
ab ab 


(9) 


Every term' at right is known here. 

On the other hand, writing V n = V • n directly in terms of the 
expressions ( 8 ) and ( 2 ), we arrive at the important relation 


ab 
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basic for what follows. 
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( 10 ) 


IV. PRINCIPLES OF WIND VECTOR DETERMINATION - 
THREE -BEAM GEOMETRY 


It will again be assumed in this section that the wind vector is 
constant throughout the observed atmosphere. 

Let us introduce a third (c-) beam identified by a point P 3 and the 
unit direction vector 


Z = 7ii + 7ai + 73 k. 
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It must not intersect with the first two beams, so that the conditions 
apply 



x 3 - x 2 

y 3 - y 2 
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( 11 ) 


Correlation of the b- and c-records and of the c- and a-records will 
result in two covariance maxima at t = t£ c and t = T* a , respectively. 

For the unknown components V x , V 2 , V 3 , three linear equation can 
than be set up according to pattern (10): 


A V 

T - ^ 2 . 
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(ap/) Vi = 7 X l x + a x L 2 + P 1 l 3 > 

(DIP/) V 2 = 72^1 + + 32^3 ) • 

( a 37) ^3 = 73 L 1 + a 3^2 + 33 L 3 > 


( 12 ) 


(13) 
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The determinant 



<*i 

a 2 

°b 

(aP7> = 

Pi 

P 2 

P3 


7 1 

72 

73 


( 14 ) 


must not be zero. None of the beam directions a, j3, 2. allowed to be 
parallel to the plane of the other two if crossed. In particular, 
parallelism of any beam pair (e.g., a ^ = p^, i = 1,2,3) is excluded; 
it would cause the determinant to vanish. 

By the system (13) three beams suffice to measure winds; however, 
not all conceivable winds. If, e.g., V happened to be parallel to the 
pair of parallel planes established by two beams, one of the right sides 
in equations (12) would be identically zero and thus remove the equation 
from the system. If V were parallel to one single beam, two equations 
would become meaningless for the same reason. Such velocities there- 
fore cannot be observed by the lineal beams as used. Although one can 
show that more than just the above few directions are not accessible 
through a set of three (single beam) detectors alone, such a set is 
being drawn on throughout the report as paradigmatic and basic. 

After solving for V l9 V 2 , V 3 the points P*, and the correspond- 
ing pairs P 3 , P* (on b and c) , P*, P* (on c and a) can be determined 
(Figure 4). In doing so, one is led to use the unit vector 
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Note that one may write V±_ for v i? since the common factor n/v^ + v| + V| 
cancels . 

Like these position vectors the lengths of the travel paths are 
needed for later application; they can be written down according to 
expression (4) where, since m = ± v, absolute marks will be added: 


p i p 2 = R ab 
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= ca 
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It is of some interest to also consider the point pairs P*, P* 
on beam a, P 3 , P '2 on b, P 5 , P 4 on c (Figure 4). Their distances 
£6 - £1 etc. appear in a mathematically compact form if the velocity 
components in expressions (15) are replaced by the solutions (13). One 
then arrives at 


£s - £l = (T ab + T bc + t m> 

. .U L 3 £ 

- li = a 


r 4 


r * 
r 5 


Ll Z 

(apy) a ' 


By conditions (1) and (11), none of the defined in the system (12) 
are ever zero. If, however, a measurement shows that a - 0, all three 
distances must have been zero, meaning that the two points on each beam 
coincide (Figure 4). An easy reasoning leads to the conclusion that in 
such a case a straight line parallel to wind direction exists containing 
the three -double points in one row and thus connecting all three of the 
beam tracts. Obviously, two of the r* will have the same sign. Evalua- 
tion of the records must have indicated, say, maxima at positive delay 
times for the (a,b)- and (b , c) -covariances , and a maximum at a (fitting) 
negative delay time for the (c,a)-covariance. 


V. LOCALLY VARYING WINDS 


We now turn to an atmosphere in which the wind vector is allowed to 
vary with location but not with time, or only insignificantly so. 
(Remarks to the latter point have been made in section III.) In other 
words, for the duration of the observation the wind motion is supposed 
to be reasonably steady in the aerodynamic sense of the word. 

With varying winds a peculiar difficulty arises in that the 
observed covariance peaks may have been caused by different winds 
blowing in different regions of the atmosphere. The system (12) must 
not be solved in these circumstances, since the times t£ c , T* a 

might be associated with three different eddy motions, V 1 from a to b, 

V" from b to c, V m from c to a. The system would break up into mutually 
independent equations. 
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To escape deception as far as possible one should keep the beams 
running close together in a narrow space region one wishes to explore. 

The wind can then be taken as constant there, and the formulism of the 
preceding section applies. Furthermore, the beams should be guided so 
as to enable a large number of unknown eddy trains to connect all three 
within that region. This furnishes us with an equally large number of 
meaningful travel time triads. If one of them is actually observed, the 
actual presence of the pertaining wind in the region is strongly indicated. 
To be sure, there might be connecting paths outside it, associated with 
eddy courses in more or less distant regions. But there will be few of 
them (Figure 3), and there is a slim chance only that the corresponding 
triads will fall within the compass of those established as meaningful, 
leading to the acceptance of a spurious wind. For added confidence, one 
may arrange for a fourth beam to be sent through the same narrow space 
volume; if the solutions of the four systems (12) that emerge with the 
possible three-beam combinations yield, by and large, four identical wind 
vectors, the likelihood of having observed a true wind can be said to 
approach certainty. 

It is seen that the problem comes up of finding serviceable beam 
configuration which can respond to a relatively large variety of con- 
stant winds that may exist in a given small space volume. It will be 
dealt with here under the simplifying assumption that the natural wind 
can be taken as horizontal. For such winds, direction more than 
strength turns out to be the quantity limiting the scope of trustworthy 
measurement. A number of trials has led to the conclusion that the 
azimuthal range should not be taken larger than 90°, or, positively 
speaking, that with three detectors horizontal winds blowing into or 
out of a quarter of the rose can be observed. This range could be 
slightly extended at the expense of mathematical simplicity (and per- 
haps physical accuracy). To cover the entire rose is not feasible. One 
cannot connect three lines by horizontal parallels which, whatever 
direction they may have, are always at about the same height. 


VI. MEASURING HORIZONTAL WINDS 


With detectors P x , P 2 , P 3 at ground level and beams pointing upward, 
z x = z 2 = 23 = 0, and Ot z >0, p 3 > 0, y 3 > 0. For ease of writing and 
greater perspicuity we introduce the abbreviations 
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7 3 
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k = 1,2. 


(17) 
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The determinants (1) and (11) can then be written as 


A ab " a ^ 3 
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= 73^3^3 


The three determinants A^ newly defined here appear in the system (13) 
for the V ± which goes into 


-iPfir) v 
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S oft ) 
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ca 




If, with a given beam configuration (yet to be developed) experi- 
mentation yields three values t*. , if , T* such that the ratio 

ab be’ ca 


v 3 


V, 


n/V? + 


+ V 


2 
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turns out to be very small, the simplification v 3 = 0 seems warranted; we 
are dealing with an essentially horizontal wind. 


On condition that V 3 = 0 the third line of the system (19) can be 

used to eliminate, e.g., the ratio A 3 /t*. . V 2 and V 2 are then given by 

at) 


^ Vx - - bj ^ + 






( 20 ) 


It will be noted here that the determinant 


Ai( c i - b i) 42 (a 1 - b i) 

A 2 (c 2 - b 2 ) A2(a 2 - b 2 ) 


A 2 A2 


( sm ) 

013&373 


is different from zero by earlier hypotheses which exclude beam inter- 
section and parallelism. Equations (20), for the moment considered as 
a linear system for 1 /t*^ an< 3 I/t^, can therefore t> e solved; the ratio 
of the travel times emerges as 


= Ik£ = ^ ( a g - b g) - q( a i - b i) 
^ - T* b A 2 q(c 2 - b 2 ) - (c 2 - b 2 ) 

where the quantity 


( 21 ) 


q 



tan cp 


( 22 ) 


defines the azimuthal angle, cp, as counted from the positive x-axis. 
The derivative of Q with respect to q shows that the function Q(q) is 
monotonic* for all q, that is, forever increasing or decreasing, with 


* 


It represents a hyperbola with a horizontal and a vertical asymptote. 
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an infinite discontinuity at q = — 1 * This particular value of q 

c l ’ b i 

is not permissible (zero or infinite travel times are not admitted); 

a - b P 

neither is the value q" - — *■ ^ , at which the numerator of Q(q) is 

a i ” a i 

zero. From physical necessity alone then certain q-ranges around q 1 
and q lf are inaccessible; horizontal winds in certain cp-ranges cannot 
be observed with a three-detector setup. 

One may arrange beam directions such that both q f and q" are negative 


c 


c 



^<0, 



(23) 


By definitions (17), these conditions impose restrictions on the direc- 
tion cosines c&£, 7 ±. 


Is it then possible to 11 catch" at least all winds with positive 
values of q(0 ^ q ^ oo ) 5 i.e., with azimuths in 0 % cp ^ 90°, 

180° ^ cp ^ 270°? In these ranges the ratio of travel times will mono- 
tonically vary between the values 


and 


Q(0) = Q 0 

Q(“) ^Qoo 


A2 3 . 2 “ b 2 

Ax c 2 - b 2 

Ae a 2 - b x 

Ai c x - b-,_ 


(24) 


These are uniquely determined by the detector positions and beam direc- 
tions. In the latter course of the investigation the quantities Q q and 
Q^, conversely, will be prescribed and then aid in shaping the experi- 
mental configuration. They constitute parameter combinations and have 
been introduced to lessen the number of significant parameters. 

To answer the above question it is at first necessary to decide on 
the location at which the measurement is to be made, especially on the 
height interval circumscribed by the positions of the points P*. The 
heights of P*, P*, P* can be disregarded as equal to those of P*, P*, 
P%, respectively. Since the z ^ are zero, the system (15) yields 
relatively simple expressions for the z- components : 
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( 25 ) 



(yg - y i ) - (x? - x i )q 
(a 2 ~ b 2 ) - (a x - b-^q 


z‘3(q) 


I y.2 - y?) - (x^ - x P )g ( 

(b 2 - c s) - (bj. - c 1 )q / ' 



(yi - y-^) - (x-, - xQq 
(c 2 - a 2 ) - (c-l - a x )q 


In a rectangular (q,z ) -system the curves z^(q) present themselves as 
hyperbolas with asymptotes parallel to the system axes (provided that 
the factors of q in the denominators are not zero). They favor our 
purpose in that they offer wide q-ranges in which the z*-values are 
nearly constant. Winds blowing at (practically) the same level from a 
rather large compass of directions have thus a chance to be detected. 
This decisive advantage is lost with the straight lines that would 
appear if the q-coef f icients in the denominators were set zero; this 
would also entail that (agy) = 0.* Let us agree then to impose the 
further cond it ions 


(a i - b x ) t 0, (b x - c x ) t 0, (c x - a x ) £ 0 


on the beams’ direction cosines. 
While it is true that 


:*(q") = ± <*, 


(q') = ± 


these discontinuities , corresponding to the vertical asymptotic lines 
q = q" and q = q* , are placed outside the interval to be investigated 
(0 ^ q ^ oo ) , both q” and q’ being negative. To ensure the same feature 
regarding the third hyperbola, z*(q), we will require that 

c^ _ a~ 

_£ 2 < o. 

C 1 ” a l 


■jV 

To avoid this predicament, only one of the coefficients in fact could be 
taken as zero. One can show that all three vanish if two do so. 
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The three equilateral hyperbolas will then, in 0 g q g oo run more or 
less parallel to their horizontal asymptotes. 

For a more detailed description, the coordinates of the three 
detectors must be taken into account, as appearing in the numerators 
of the functions z'^(q). Let us place P x at the origin, P 2 and P 3 in 
the upper right and the lower left quadrants, respectively (Figure 5). 
This arrangement has been found advantageous for measuring horizontal 



Figure 5. Detector arrangement and coordinates. 


winds with azimuthal angles in the remaining two quadrants. The four 
positive position parameters x 0 , y 0 , A, (j. are introduced for practical 
reasons and are left indeterminate at present. With their aid the set 
(24) can be rearranged to give 


zl(q) = 


z 3 (q) = 


z*(q) 


m y + q 
o 


(a 2 “ b 2 ) ~ ( a i ~ ^jl)^ 
(li + 1 )y o + (A + l)x Q q 

(b 2 - c 2 ) - (b x - c x )q 


V + X a 




a 2 ) - ( c i - a x )q 


(25) 
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These heights all ought to be near a preselected height, z = h > 0, near 
which the constant wind vector is being sought. This is especially true 
at the proposed ends, q = 0 and q = oo, of the q- interval. With positive 
numbers D x , E^ not far from unity we may therefore put 



hDi, 


(n + 1 )y f , 

- - bTHsf = 



(q = o) 



hEi, 


(A + 1)x q 

b l " c i 


hE 3 , 



hE s (q 


(26) 


As an immediate consequence 


( a 2 “ b 2 ) > 0, (b 2 - c 2 ) < 0, 

( a i " b j) < 0 , (b - c j) > 0 , 



(27) 


These conditions summarize and specify all restrictions placed so far 
on the direction cosines. They still leave much room for free choice. 
But, if they are satisfied, the values of z*, z*, z*, for all wind direc- 
tions 0 ^ q ^ oo, vary monotonically* between hD x and hEj, hD 3 and hE 3 , 
hDs and hEs, thus insuring that one is measuring effectively "at" the 
height h. This, however, is not enough. 


One must also see to it that, near height h, the three beams are 
close to each other; i.e., that the travel paths R^b> ^bc> Rea are 
reasonably short line segments. Evaluation of expressions (16) in terms 
of horizontal winds and the quantities introduced in the meantime yields 


ny r 


R* = sll + q 2 
ab jov E_i 

A D x 


Pi - Ex 




D- 


«b 


f + u - y ° 


u + 1 E3 . 

tn v D? + q 


~ 


D- 


R* = n/T + q 2 
ca n 


y 0 




D 5 z £ 5 


(28) 


*This can be seen by consulting the differential quotients dz^/dq. It 
is also evident from the course of a hyperbola near its horizontal 
asymptote. 
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where v stands for y 0 /x 0 : 


v 




( 29 ) 


These results indicate that, even though the D and E are all to be chosen 
close to unity, one must select them such that 


2l = . a .3, h. A i 

E^ A 3.2 ~ bg 


2a 

e 3 


u + 1 
V + 1 



1 >. 


(30) 


2s 

Es 


*“ V 


~ + l 

C 2 “ 3-2 


J 


More conditions for the relationship of direction cosines and detector 
locations will appear in section VII. 

At this station it will merely be pointed out that the distances 
(28) attain minima for certain values q = q min : 


r, 

ab 


R* 

R b< 


ca 


A 


with 


> 

= X 

min ° /l + q^. n 

Ei 

c . 

4 min )j,v E ! 


A + 1 

d 3 - e 3 

with 

A + 1 

I 2a 

- X 

min 0 vl + q 2 . 

min 

e 3 

"min (i + l 

V E 3 


1 

D s - E 5 

min ° n/ 1 + q^ in 

e 5 


1 £= 

with q . = — — 

min v E 5 


(31) 
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While at q = 0 and q = oo the points P /v must not be too far apart to avoid 
physically prohibited long travel times, the minimum distances must not 
be too short to prevent experimental inaccuracies from falsifying the 
outcome.* 


The core of the space volume traversed by closely bundled beams is 
a triangular prism with (non-parallel) bases, P* P* P* and P* P* P* on 
figure 4. Their shape is not critical; no further restrictions on 
direction cosines issue from position requirements for the corners. 

The latter's heights (z-coordinates) are approximately equal. Computa- 
tion (with v 3 = 0) of the x- and y- coordinates as appearing in the 
formulas (15) reveals that they all can be written in terms of z*, z*, 


X* = 
X 1 

a i z i 

(q) 

A- 

a 2 z*(q) 

X* = 
X 2 

-Ax o 

+ b 1 z , £(q) 

vl- 

^y 0 

+ b 2 z*(q) 

K 
u * 

ii 

-Ax 0 

+ b-jZ^Cq) 

y 3 " 

ny 0 

+ b 2 z*(q) 

x * = 
X 4 

x 0 + 

c i4(q) 

y * ■ 

-y Q 

+ c 2 z*(q) 

x * = 
X 5 

x o + 

Cl z*(q) 

•a. 

y 5 = 

■y 0 

+ 

* 

X 6 = 

a i4 

(q) 

y| - 

a 2 z 

5(q) 


Although these coordinates depend on q (i.e., on wind direction), they 
ordinarily would not vary very much, as the z‘£(q) remain close to the 
fixed height h. An exception occurs if one wishes to examine a small 
space volume at a large lateral distance from the origin. The beams 
then must be sharply inclined towards the ground causing some of the 
l a kl> l b kl> | c k I to he large, perhaps much larger than unity. In such 
cases the lateral location of the space volume cannot be adequately 
known beforehand, because it will markedly shift with the still 
undetermined wind direction prevailing at height h. By contrast, if, 
e.g., the a-beam is sent up vertically (a x = a 2 = 0) , one is sure to 
measure more or less overhead, the corners P* and of the prism being 
located on the z-axis. The other four corners remain dependent on q, 
however. 


If from an independent source the wind is known as horizontal, these 
requirements could be relaxed regarding R“ a , since is not needed 
for the determination of V z and V^. 
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In concluding this section the system (20) will be put into a form 
more suitable for the argumentation in the next one. By repeatedly 
employing expressions (24) and (26) one finally arrives at the 
representation 


FV X = 


FV 2 = 





(33) 


where 


F = 


( a l “ b x)( c 2 ” b 2) " ( a 2 “ b 2) ( c 1 " b l) 


1 - -£ 


c i - b i ^ 0 ( a 2 - b 2 ) + My Q ( a i - b i) 


_i 

Ax D 


% 

Q, 


X- 1 


C = 


= ^ ■■ b ' 2 = f±4 v ^ > 0 

b x -c i A + 1 D3 




(34) 


The ten parameters x Q , A, p, v, a^, b^, c^ are here reduced to four: 
F> C, Qoo, Q 0 . 

Looking at F it is clearly necessary that 


a i “ b i c 2 _ b 2 

Q c C 1 “ b l a 2 " b 2 


(35) 


By requirements (27) this ratio is positive, so that and Q q must have 
the same sign in addition. That Dx/E x should be different from unity 
has already been deduced from the physical necessity to have the travel 
paths as different from zero; see set (30). The lines a and b must 

not intersect. 
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Truly important are the last three parameters, as will be seen 
shortly. If they can be assigned values so as to keep error trans- 
mission, height variation, and detector distances in bounds, the problem 
will have been solved. Pertinent analyses are offered in the next two 
sections. They will show the requirements as partly incompatible; a 
trade-off is unavoidable. 

The knowledge of all four parameters is necessary for computing V x 
and V 2 from the set (33). 

VII. ERROR PROPAGATION IN DETERMINING HORIZONTAL WINDS 


Freedom in the choice of the beam arrangement is further narrowed 
down by the desire to minimize the effect of observational errors on 
the wind vector results. For the investigation to follow let us agree 
on a pivotal travel time of 1 second subject to a 10 percent error in 
ascertaining it. The ensuing uncertainty of ±0.1 second will be taken 
as typical for all travel times. This excludes r f s materially smaller 
than 1 second as suspect of too great inaccuracy. The observational 
error should remain small relative to the peak times.* 

If winds up to k knots, i.e., up to (approximately) 1/2 k m/sec, 
are to be measured with some confidence, none of the minimum path lengths 
(31) will be allowed to be appreciably shorter than 1/2 km, so that not 
less than 1 second is needed for the eddies to cover the travel distance 
in between beams. This places a restriction on detector locations and 
beam directions. Another one stems from the necessity of limiting the 
lengths of the longest paths (at q =* 0 and q = oo) in order not to 
obliterate the correlation. It is bound up with the behavior of the 
turbulent flow under observation, and therefore less definite unless 
pertinent flow characteristics can reasonably well be guessed at. 

By the system (33) the magnitude of the wind vector becomes 


V 




(36) 


Other assumptions on it would of course alter the analysis of the 
present section. 
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The direction of the wind is given by 
V i Q ~ Q ro 

— = tan cp = q - q q~Tq- ( 3 7) 

No decision can be made here on the (equal) signs of and Q q . Whether 
they are both positive or both negative, the condition 00 


0 ^ q ^ oo 

is always met, the constant C being positive. 

There are practical advantages in using a beam configuration where 
T^ b and t§ c are of equal sign (Q >0); the mathematical notation is 
simpler, and the records handle more easily, since the travel time r* a 
will as a rule have the opposite sign. To show this, suppose that the 
beams are running such that, if at height h a horizontal wind is blowing 
with an azimuth in 0 cp ^ 90°, the registering eddies proceed from 
a to b (T* b > 0) and from b to c (t£ c > 0) (although at slightly dif- 
ferent levels close to h) . At a third level near h they should move 
from a to c in a majority of cases; but, since in Figure 4 the "positive" 
direction m was taken from c to a, m = -v, causing x* a to be negative 
(compare expression (5) for V, as it would be written for the transition 
c -> a) . 

Let us then simply require that 

0 < Q g Q £ Q . (38) 

oo o 


Is there a way of assigning favorable values to the parameters 
Qoq, Q 0 , C such that variations of ±0.1 second in the t 1 s have as little 
effect as possible on V and cp? Those values could not be expected to 
minimize the errors in all admissible wind vectors; in fact, they can 
achieve that for some directions and strengths only. However, one may 
endeavor to keep the error within certain bounds, say 10° in direction 
and 107 o in strength. 

As regards direction errors it follows from relation (37) that 


dcp 


cos 2 cp 

C 



- Q 


- Q) 


2 


dQ. 
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Writing, for short, 


T*_ = T ls 


Jc _ 


ab ‘ l5 ‘be T2 ’ 

T l dx 2 " t 2 dT x dx 2 - Qdxi 


dQ = 


T l 


(39) 


Thus, in a first approximation, the angle error engendered by the errors 
At 3. and At 2 becomes 


c(Q 0 " <U at 2 - QAti 

^ = C^(Q q - Q ) 2 + (Q - Qj 2 “ 
Since by definition (21) 

Ti T 2 ’ 


(40) 


(41) 


it is seen that, at a fixed value of Q, 

| x x | or |t 2 | decrease. This is further 
times below 1 second. 

The observational errors have been 
ing on the worst case we have therefore 


the error will increase when 
reason for not admitting travel 


limited to ± 0.1 second. Draw- 
to consider the error possibilities 


At ! = At 2 = ± 0.1 sec 
At x = -At o = + 0.1 sec 


so that 


= AT2 


C ^ 0 ' Q oo> 


Q ) 2 + (Q - 




1 + Q 

Tl 


(42) 


(43) 


With the upper sign the errors Ap are relatively small. Concentrating 
on the second possibility we shall take either t x = +1, or t 2 ~ +l x 
(the lowest values admitted). 


^Negative unit values for t^ or t 2 result in the same errors with the 
opposite sign. 
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Dealing at first with t x = +1 and writing 


api = At 2 • C (Q 0 - QjFiCQ), 

where 

Fl(Q) = c 2 (Q 0 - q) 2 °+ (Q - q m ) 2 ’ 

one is restricted to the interval 1 g Q ^ corresponding to t 2 = 1 . 

Taking secondly t 2 = +1, one has to write 


= At 2 • C(Q 0 - Q oq )F 2 (Q) 

with 

•p ro >, = Q(! + Q) „ 

2(Q) C 2 (Q d - Q) 2 + (Q - Q ro ) 2 

and the restriction i Q ^ 1, as ^ 1 here. Note that to set t 2 = 1 

in expression (43), one must introduce in it the relation (41). 


On differentiating the functions F x and F 2 with respect to their 
argument Q one finds that both can have maxima within the respective 
Q-intervals on condition that 


1 - Qoo Qoo + 3 

c2 = K 1 - q-^T 

o o 


in the first case, 


C 2 § K 2 


1 - Qco 3 Q~ + 1 

Qo - 1 + 1 


in the second case. 


(44) 


Neither the values Q* and Q 2 at which these maxima reside, nor the corre- 
sponding maximal errors are given here, because, with the relationship of 
Qoo> Qo an< 3 C 2 adopted later, the above conditions are not satisfied. In 
these circumstances the largest absolute error occurs at an interval 
terminal; one can show that in both cases this is Q = 1 . The errors 
and £cp 2 are then equal and amount to 


largest 


1. 

5 C*(Q 


C (Q 0 - Qco) 

o - + <! - V 


> 


(45) 
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At 2 being taken as ± 1/10. Notice that this is an estimate only for the 
actual error (differences have been treated as differentials). 

Up to now nothing more definite can be said about the three parameters 
than that they should so be selected that the largest angle error won* t 
exceed 10°. Further hints would be expected to issue from a discussion of 
the strength errors. 

With the aid of relation (41) logarithmic differentiation of expres- 
sion (36) for V yields 


dv 1 (Q " Q oo> ( -Q 2dT l + Q co dT 2) + C 2 (-Q + Q 0 )(Q 2 dTl - Qo d T 2 ) 

T “ t 2 (Q - Qoo ) 2 + C y (- Q + • 

If again differences are written for differentials, the relative strength 
error in terms of Ax! and Ax 2 is being estimated as 


AV 1 Q 2 At 1 [C 2 (Q o - Q) - (Q - Q ro )] + At 2 [Q co (Q - Q*) - Q q C 2 (Q o - Q)] 

V “ t 2 (Q - Q m ) 2 + C 2 (Q o - Q) 2 

(46) 


Like the cp-error it changes sign with t 2 (and Ti) and absolutely increases 
when Tj_ or r 2 (taken as positive) decrease. 

We wish to determine the worst possible errors. With this end in mind 
let us consider values of Qoo and Q 0 such that 

Qoo < 1 < Qo- < 47 > 


Then if xi has the smallest admissible (worst) value (x^. = 1), the travel 
time x 2 can vary in 


1 - T 2 - t iQ 0 “ Qq' 


(48a) 


if x 2 is left constant as unity, Tj can move in 


t 2 i 

1 " Xl ' Qoo " Qco * 


(48b) 
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It is advisable to examine first the second of the possibilities 

(42): 

/AVs Ax 2 <Q - QcoXQoo + Q 2 ) - c 2 (q o - Q)(Q o + Q 2 ) 

(v) 2 - T2 (Q - + C*(Q 0 - Q) 2 * (49) 


At Q - Q m , t 2 can be taken as unity; can be taken so at Q = Q 0 . Hence 

the errors at the Q-interval boundaries will be written as 




<— > 

^ V } 2 


Q = Q, 


AT 2 Qq + 

’ T 2 Qq - Qco 

\ 

AT2 Qc + Qq ( 

T 1 V Q 0 ‘ 


(50) 


In the last expression relation (41) has been put to use again. 

Assuming Q q » and |at 2 | = 0.1, the worst terminal errors (at 
t 2 = 1 and' * 1, respectively) are of the order of 10 percent. Since 
they are of opposite sign, the error zero occurs at least once in 

^ Q ^ Q 0 (the observation errors cancelling each other). The value 
Q = 1 with which both and r 2 can be unity belongs to both the intervals 
(48a) and (48b) (standing at the lower end). While nothing much can be 
done to reduce the error at the upper end, we can see to it that at Q = 1 
it becomes zero. Expression (49) then yields the relation: 



It simplifies further if we require (as seems best to do) that the longest 
travel paths associated with Tj_ = 1 and t 2 = 1 should be equal. These are 
covered in Q q and 1/Q^ sec., respectively, as can be seen by the ranges 
(48). Thus, 


Q 0 - 1/Q„ (52) 

and 

c = Qco. (53) 
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As a consequence of these stipulations the worst relative error (49) 
may now be written as 


A QJQ 3 - i> - Q e + Q "l 

l »'s t 2 Q to (Q^ + 1) - 2Q 


where 


- 


Q + Q 

O x oo 



T 2 = 1 in Qoo - Q - 1 

t 2 = Qt x = Q in 1 g Q g Q q 


(54) 


Clearly, the relative error is 0 at Q = 1. No other zero error (AV/V) 
exists, since the equation 2 


Q m (Q 2 + Q + 1) - Q = 0 


has no real roots. 


Relation (37) now attains the form 


q = 



00 


00 


and shows that q = 1 when Q = 1. That is, we have decided that the strength 
error zero existing when At x = -At^ should occur for a wind with azimuth 
45° (or 225°). 

In the event that At x = At 2 (first of the possibilities (42)) one 
finds that 


M , at s -VQ 3 + 1) + Q 2 + Q 
1 ' t 2 <^(Q y + 1) - 2Q 
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Q = Qco 


( 56 ) 


(-) 
v v'l 



l 


Q=Q 0 


2 "X 

At 2 Q 0 ~ Qoo 
' ~ Q 0 - Q ro 

^ * 

At 2 q M ~ Qq 
T 1 Q 0 (Q o - Qoo) „ 


In formula (55) the relationships (52) and (53) are incorporated. The 
worst terminal errors again are about 10 percent with Q q » Q^, both of 
the same sign this time. One can show that there is no zero of the 
function (AV/V^ in Q M g Q s Q 0 ; however, absolute errors less than 
those at the terminals will prevail in the interior of the Q- interval. 

Turning to the angle error we notice that the stipulations (52) and 
(53) can be combined to 


C = 



so that, with < 1 , 

Ki Q o " 1 Q<x + 3 Qoo + 3 

C 2 = 1 - Qco Q 0 + 3 = 1 + 3 Q m > lj 

rr 3 Q~+ 1 

= < i 

c 2 3 + < A * 

Conditions (44) are violated; the error expression (45) applies and can 
now be written as 


il = JL 1 + Q ” 

'^1 largest 10 1 - * 

If the largest deviation permitted is 10° ~ 0.175, one obtains an upper 
bound for as well as a corresponding lower bound for Q q : 


Qco - n 


a y = 3.67. 


For deciding on a definite value for Q q we recall that the ratio of 
travel times Q for a (constant) wind is equal to the ratio of the eddy 
travel paths. Thus, with cp = 0, 


X 2 


q=0 


be 

* 

ab 


q=0 


Since 


ab 


q=0 


is already longer than the minimum path, one would not care to select Q q 
too large. Table I below (slide rule computation) is based on 

Qo = 5 - 

If a turbulence pattern moving parallel to x-axis needs 1 second to cover 
the distance between the a- and b -beams , it requires 5 seconds for the 
distance between the b- and c-beams and therfore is assumed here to essen- 
tially preserve its characteristics during any 5 seconds of of observation 
time in order to insure a well-developed covariance peak. With a speed of 
20 m/sec this corresponds to a 100 m stretch. 

It is important to keep in mind that and Q q (corresponding to 
cp = 90° and cp = 0°) are the limits of serviceable peak time ratios. 
Observed ratios 





and 



> 5 


ought to be discarded since there is suspicion that they might not have 
been produced by a single wind present near height h. 

For computing Table I the now particularly simple and symmetric 
equations (33) are available: 



^2 
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The indices 1 and 2 are interchangeable without altering the form of the 
equations. It suffices therefore to deal with the half-range UQ gQ 
only; the same absolute errors prevail in ^ Q g 1 . The value of 5F ° 
is inconsequential for our purpose which is concerned with ratios only: 



The primed quantities assume that the peak times as determined by the 
records deviate by as much as ±0.1 second from the true times, and t 2 
While these are already the largest deviations considered likely to occur 
the error situation is still worsened by adopting the smallest values, 

Ti = 1, Ti = 1.1, as can be done in the above half-range. These are 
implicit in Table I which explicitly shows the and T^-ranges only. 

TABLE I 

Worst Expected Errors (Q q = 5) 


^2 

*2 

- 1) X 100% 

0- 

1 

*0- 

I 

— 

(~ - 1) X 100% 

cp' - cp 

5.0 

4.9 

-9.5 

1 °28' 

5.1 

-9.5 

0°56' 

4.5 

4.4 

-9.6 

1 °41' 

4.6 

-9.5 

1°1' 

4.0 

3.9 

-9.4 

1 °57' 

4.1 

-9.4 

1°7' 

3.5 

3.4 

-9.3 

2°19' 

3. 6 

-9.2 

1°10' 

3.0 

2.9 

-9.2 

2°51' 

3.1 

-9.3 

1 °21' 

2.5 

2.4 

-9.0 

3°35' 

2.6 

-9.1 

1 °31 ' 

2.0 

1.9 

-7.8 

4°46' 

2.1 

-8.8 

1°28' 

1.5 

1.4 

-5.3 

6°30' 

1.6 

-8.5 

1°13' 

1.1 

1.0 

0 

8 °10' 

1.2 

-8.6 

0°13' 
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The above errors are found little affected by the choice of Q 0 as 
long as Q ^4. Attempts to reduce them by substituting for the condi- 
tions (51) and (52) other relations judged promising* have met with no 
decisive success; improvement in one place had as a rule to be paid for 
by worsening elsewhere. 

The errors shrink when both travel times are larger than pivotal. 
This reflects the relative smallness of the constant observation error 
in such cases and could not be expected to occur if the latter, for 
example, were to be taken as increasing proportionately to observed 
peak time. 


VIII. EXPERIMENTAL ARRANGEMENT 


The (relatively large) value selected above for Q 0 cannot be main- 
tained. It must be relaxed in order 

(1) to keep the D^, Ej_ close to unity, 

(2) to avoid too sprawling a detector configuration, and 

(3) to secure reasonable travel path lengths. 

The peak errors then will go up. We shall strive to hold them below 20° 
in direction and 15 percent in strength. 

It will be seen shortly that conditions (51) and (52) need not be 
altered . 

To inquire into these points we will have to make use of relation- 
ships that exist among the position parameters A, p, v, the bounds Q^, 

Q , and the D^, E^. 


Vc 

Such conditions were derived from the structure of the error approxima- 
tions (43) and (46), and similar ones. Among others the case 1 < < Q q 

was s tud ied . 
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The first line in the set (26) may be written as 


a 2 “ b 2 — 


i^c 

hD- 


b 2 - c 2 - 


(n + i)y c 

hD^ 


C 2 “ 3 . 2 — 


hDc 


Since the left sides add up to zero, 


( 57 ) 



D 


D 3 - D 


3 

1 


n + 1 


3 

D 5 



In an analogous manner, the second line yields 


(58) 



e 5 - e 3 

E 3 - E x ’ 


A + 1 


h 

E 5 



(59) 


Because A and p are defined as positive, D 3 , D 5 and E x , E 3 , E 5 must 

be in either ascending or descending sequence. 

Again using the first line, we may write the first expression (24) as 


Qo 

Both the lines 
into 


u £3 Ag 

(26) help to transform 


the definitions (18) for A^ and Ai 


x 0 y 0 1 1 

(A+iXu + D^- (---) 


- 


Ap 


x 0 y 0 

h 
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s o tha t 


Q o 


A + 1 E T P 3 ~ E 3 

A E — D E 3 


Analogous transformations yield 


(60) 


Qco - - 


A £5 ^2 = j£ + 1 °3 ~ E 3 Pi 
A + 1 E x A x [i D 3 E x - D x 


(61) 


Solving equations (60) and (61) for D^E! and D 3 /E 3 one obtains 


1 + ^ 
r> A Q 

£l _ 

Ei ~ ! + M_+. I -L 

H Qco 

A „ 


(62) 


£2 

E 3 


1 + 


A + 1 


1 + 


— q 

(1 + 1 w °° J 


By expressions (31), these quotients exercise the dominating influence on 
the minimum travel paths and R£ c * Which one of these is the shorter 

one cannot be seen at a glance. To clarify this question we note at first 
that by relations (34) and (62), 


Ha 

e 3 


u + 1 v_ 
A + 1 C 


Hi. 

A + 1 

H Q “ _ 

Q 

j I 00 V 

Ei 

A 

H + 1 Qo E 3 

AQo c 

With the aid of the results (62) 

and (63), 






\c 


1 C 2 ^ 2 + Q 2 

f o 00 


R * 

ab 

min _ j 
min 

C 2 + 1 ’ 



(63) 
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This root can have positive values smaller than, larger than, or equal 
to unity. It would still have to be compared to (more involved) analogous 
ratios formed with 



min. 


The virtually infinite number of possible minimum lengths is drastically 
reduced if one limits the competing cases by stipulating the paths 




and 


min 


R 


ab 


min 


as equal, i.e., by requiring that 



Surprisingly, this condition* is the same as was found earlier in demand- 
ing that the zero strength error should occur at Q = 1. There is no 
reason why the second of the former conditions should not be adopted 
as well: 


Q 0 Q°° 1 > 

so that the simple relation 


C 


= Q 


00 


( 52 ) 


( 53 ) 


again applies. 

To keep the minimum path lengths different from zero none of the 
ratios Di/E^ must be equal to unity. On the other hand, the beams, as 
is required of them, will be running at approximately the same height in 


It means that the equal minimum distances are found between the heights 
z x (q) and z 3 (q) defined by q = and q = Q 0 (both of these values are 
in 0 S q g a.) . 
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parts if the and Ei are allowed to range between, say, 0.9 and 1.1 
only. Their ratios then move between 0.818 and 1.222 with exclusion of 1. 

For a more detailed investigation take 



to start with. The first of the expressions (62) coupled with condition 
(52) then yields 

^<■ 4 ^- (64) 

With \i > 0 (as was agreed upon) the right side must be positive: 

(65) 

From Figure 5, 

I . _ 

A x 2 5 

so that Q = 5 would require an uncomfortably large value for this ratio 
of P 3 , P 2 -coordinates . 

In the same way, when limiting D^Ei to an upper bound of 1.2, we 
obtain the further condition 


1 _ 


>! 


A + 1 

. A 



1 . 


( 66 ) 


If we tentatively introduce the values Q q = 2 (Q M = 1/2) and 1/A = 5 
compatible with the unequality (65), the quantity p = according 

to inequalities (64) and (66), should be taken out of the range 2 < p. < 6. 
Let us select p = 4. 
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Before proceeding further, it is advisable to probe into the (approxi- 
mate) error peaks to be expected when using the much lowered value of Q 0 . 

Regarding the angle error," 


largest 


1 + Q 

_1 OO 

10 1 - Q 

OO 


0.3. 


This corresponds to 17°11' . 

It suffices to compute the worst strength errors at Q = Q 0 . With 
|At 2 | = 0.1 and = 1 expressions (50) and (56) give 


K V*2 


Q=Q r 


- ^ (~ - 15 %)** 


A 
' v ' l 


Q-Q, 


_7_ 

60 


(~ -11.7%). 


This appears bearable. A table containing more accurate error 
results is supplied at the end of the section. 

Based on Q 0 = 2, A = 1/5, p = 4, 



1.14 



min 


2x 

o 

35 n/7 


0.0256x = R* 

O DC 


min 


/V 

Conditions (44) are not satisfied, as they were not with Q q = 5. 

' To comply with the later Table II, we have taken here Aix =0.1, 
At 2 = -0.1. 
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The factor of x 0 would increase (as one might desire), were the parameter 
[X to be taken as 5, for instance; but then the third distance. 



min 


would become markedly shorter. For calculating it we have to use the 
quantity 

y o A + 1 4 

V "* C — if ■ 35 • 

where C = = 1/2; D 3 /E 3 is found from the second of the expressions (62). 

With the smallest admitted value, 9/11 = 0.818, for D 5 /E 5 one obtains 



min 


0.0252x o , 


which value increases to that of 


* 

ab 


min 


on relaxing D 5 /E 5 to 0.816. 


A wind of 20 m/sec (~ 40 knots) then covers a minimum distance in 
one second (smallest travel time permitted) if 


x = 780 m. 
o 


(Lesser winds require proportionately less.) Hence, the detector coordin 
ates must be chosen as 


Pi: 

xi ■ 0, 

0 

II 

H 

P 2 : 

O 

& 

1 

II 

c\] 

X 

= -156 m, y 2 = |a.vx Q = 356.7 m. 

P3: 

x 3 = X Q = 

780 m, y 3 = - vx q = -89.1 m. 
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These figures represent the most compact detector configuration that can 
be obtained with the values chosen for Q 0 , A and y., if winds up to 40 knots 
are to be measured reliably. 

A check into the longest paths (at q = 0 and q * oo) with the aid of 
expressions (28) reveals that all beam distances remain below 44.6 m, 
except the distance R* a which drops from 143.5 m at cp = 0 (q = 0) to 
47.2 m at cp = 7° (q = 0.123) and continues down to the minimum. The 
relatively long c-a-paths may sometimes cause difficulties in covariance 
peak definition, but on the whole the beams appear running sufficiently 
close in the region proposed for measurement. 

If the actual eddies are suspected to be very sizable, one might 
wish to increase the minimum distance to avoid spurious or untrustworthy 
correlation maxima. The detectors must then be located farther apart 
than above. 

It will be noted that the observation height h does not directly 
enter into the determination of the detector coordinates.* However, 
according to set (26), it affects the calculation of the beam direc- 
tions which also depends on the quantities D-^ and Ej_. 

The small value assigned to the ratio D 5 /E 5 suggests using the 
lowest possible figure for D 5 . If we put 

D 5 = 0.9, 

E 5 = 1.102.** 

The rest of the D-j_ and E^ is found in a more cumbersome way. By com- 
bining expressions (58) through (61), one obtains the expressions 


H 

A 


Ei Es 

Ei D 5 


3 1 

Es ~ E3 

P3 I Pi 

d 5 - d 3 



1 + 


1 + 


E 5 - E 3 

- P i 

d 5 - d 3 


An indirect influence will be recognized later. 


This figure, slightly larger than the admitted value 1.1, arises 
through adopting the ratio 0.816 instead of 0.818. 
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After applying the second of relations (63) and solving for 


5 = 


i _ 

D 3 - Di Q o 


Dc 


'3 C 5 

--T - 1 

v 6 5 


e = 


E 3 - E; 
E 5 - E- 


O 


So C ^ _ 

Qoo V E 5 C 


D c 

- — - 1 

v E 5 


one uses expressions (58) and (59) to obtain 


Di = mD 5 & 


E “ AE 3 £ 


d 3 - (m- + i) D 5 g + ^ e 3 - (A + l) e 5 e + ^ 


With the parameter values chosen so far, 

* 3_ _ 10.71 

0 " 10.28 ’ 2.57 


Dx = 


2.7 

2.57 


1.051 


= 1-89 

1 2.056 


0.92 


D 3 = 


13.5 

13.28 


1.017 


E 3 


5.67 

5.312 


1.068 


D 5 = 0.9 


E 5 = 1.102. 


The D's and E*s are in descending and ascending sequence, respectively. 
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These results imply that we have to require wind constancy from 0.9 h 
to 1.05 h at cp = 0, and from 0.92 h to 1.1 h at cp - 90°. In between the 
requirements are less. The three hyperbolas z*(q) approach each other as cp 
increases from cp = 0 on and actually intersect at one value* of q = tan cp; 
for it, P'£, P3, P^ lie on the same horizontal line, rendering the required 
layer thickness zero. Further increase of cp widens the traces up to the 
distance of 0.18 h at cp = 90°. 

If one intends to measure at a large height, the vertical layer near 
the cp- terminals may be judged too extended for the winds to be constant in. 
In such a predicament, one would calculate the stretch, q x ^ q ^ q 2 of 
the zi-curves in which they can be considered sufficiently close. This 
would exclude the measurement of winds blowing in a certain range of azi- 
muths near cp = 0 and cp = 90°. Mathematically, new Q-terminals, Qi > Q , 

Q 2 < Q 0 , obeying the relation 


Qco + C< ll,2 Q 0 

Ql ' 2 " 1 + Cq 1} £ 


Qco + qi >2 
+ Qoo qi,2 


would be introduced, causing a curtailing of the admissible T-ratios and 
therefore of measurable wind directions. If this is undesirable, one 
would have' to keep the and Ei closer to unity than before with the 
consequence of a more sprawling detector layout, especially since, in 
addition, higher up one should have to allow for strong winds. 

If we plan to use the detector locations as given above for measure- 
ment at h = 100 m, wind vector constancy is required at most from 8 m 
below to 10 m above h (when cp = 90°), and from 10 m below to 5 m above 
(when cp = 0°). Assuming this as warranted we can proceed to find the 
beam directions. 

The system (57) for computing a 2 , b 2 , c 2 is underdetermined, since 
there is a linear relationship of the left sides. As a consequence, one 
equation must be left out. There is still freedom in choosing one of the 
unknowns, say a 2 = OL^jOLj,. Thus 


fe _ 

a 2 

4 

X 89 

CL 2 

P3 = 

013 

100 

X 1.05 

0^3 

7£ _ 

a 2 


89 


73 " 

a 3 

100 

X 0.9 

a 3 




c 2 = 


>v 1 

This is q = 0.275, corresponding to z* = 1.035 h. 
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In the same way, the second line (26) gives 


Pi 

P 3 


Ax 


a x 


b ^~ = ai + iil " +1 - 70 


7 1 x o «i 

ex = = a i ■ rr - = ~ - 7.08. 


7 3 


hEc 


a 3 


Together with the orthogonality conditions 
expressions determine (3 and once a is chosen. Take, for example, 


a? ■ 1 , etc. ) , these 


a 3 = 7 (73.4°), 


CC 2 = y (64.6°), a± = y (31°). 


Then, 


P3 = 536 (7S. 8 °), P 2 = 


-HI (111 - 5 °>» pi-HS (24 ' 3 °> 


73 = 08 ( 78 * 2 °)> 7 2 “ of (3 9 .4°) ; 


4.08 /->/£. io\ 
” ' 08 (146 ' 7 } 


The figures in parentheses give the angles made by the three beams with 
the positive coordinate axes. Those with the z-axis are rather large; 
the beams have relatively low inclinations towards the ground plane. 

When they approach each other in the region of interest, many of their 
horizontal connections can therefore run approximately at the same 
height (those with directions in 0 ^ q ^ oo) . 

Flatness of at least two beam courses is required at any height one 
wishes to explore. If it is large, the detectors will have to be placed 
widely apart, or else one must be content to monitor a smaller sector of 
the rose. (It will be remembered that covariance peaks could originate 
through unrelated winds and that, in order to minimize that danger, one 
has to force the beams into close neighborhood near h so that the atmos- 
pheric motion there can be considered uniform.) 

Actual analysis errors encountered with Q q = 2 are given in Table II 
which again is based on t x = 1 , T x = 1.1 and is computed only for 
1 ^ Q ^ Q q = 2. It is assumed as before that the worst observational 
errors in t 2 are ±0.1 second. Observed travel time ratios t 2 /t x > 2 
and < 1/2 ought to be discarded on suspicion they might not be related 
to a single wind blowing near the height h = 100 m. 
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TABLE II 


Worst Expected Errors with Q q = 2 



i 

- 1) X 100% 

cp' - 9 

t 2 

- 1) X 100% 

9 ' - cp 

2 

1.9 

-13.3 

6.4° 

2.1 

-10.4 

1.9° 

1.6 

1.5 

- 12.1 

9.9° 

1.7 

- 10.0 

2 . 0 ° 

1.3 

1.2 

- 7.6 

14.0° 

1.4 

- 9.0 

1 . 6 ° 

1.1 

1.0 

0 

16.2° 

1.2 

- 8.9 

0.7° 


These errors reduce with slower winds, i.e., with larger values of 
the t* s. A wind half the strength of that in Table II would (with = 2, 
t 2 = 4) be subject to a largest error (incurred with At x =0.1, At 2 = -0.1) 
of -7% in strength and 3° in direction, as compared to -13.3% and 6.4°, 
respectively. In like circumstances the large figure Ap = 16.2° on the 
last line would reduce to Ap = 8.2°. 

In general, if one decides not to use observed travel times below, 
say, 2 seconds (instead of 1 second), the requirements for the experi- 
mental arrangement will be less stringent under the assumption made in 
the present report of an observational error |at| ^ 0.1 second. The 
layout can be modified accordingly. The one described here is but one 
example; others can be constructed following the same guidelines with 
shifted accents or altered assumptions. A simple case in question is to 
select a different free beam direction. 


IX . CON CLUD ING REMARKS 


Of the many parameters entering the problem only two remained 
arbitrary at the end, both direction cosines. Perhaps use could be made 
of the discretion thus given to us to forge a link with a fourth beam 
destined either to widen the range of azimuths measurable at one sitting 
or to explore winds at a different height simultaneously (four-detector 
arrangement). Similar aims could be pursued by utilizing existing multiple 
detectors that receive light from several directions (residing in one plane, 
however, with present versions). 
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If general (non-horizontal) winds are to be detected, the mathematics 
are considerably more involved. The number of height functions z* will 
increase from three to six. They will depend on two variables (q = v 2 /v x 
and p = v 3/ v x ) and thus describe six surfaces instead of the three rela- 
tively simple hyperbolic lines. A rather large effort undoubtedly will 
have to be spent in solving the problem. It might be preferable to at 
first deal with winds, outside of storm clouds for example, that are 
predominantly vertical rather than horizontal. 

Detectors mounted on a plane or satellite flying at constant veloc- 
ity present different complications. Mathematically, two Galilei systems 
moving relative to each other are necessary for the description; physically, 
the limited choice of detector positions hampers the freedom in adjusting 
the parameters x G> A, (_t, v to meet given requirements. Some of them 
probably will have to be relaxed. Which to select and how far to yield 
will have to be discussed in the course of later work. 
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